ANALYSIS OF A MODEL ARISING FROM INVASION BY 
PRECURSOR AND DIFFERENTIATED CELLS 



XIAOJIE HOU 



Abstract. We study the wave solutions for a degenerated reaction diffusion 
system arising from the invasion of cells. We show that there exists a fam- 
ily of waves for the wave speed larger than or equals a certain number, and 
below which there is no monotonic wave solutions. We also investigate the 
monotonicity, uniqueness and asymptotics of the waves. 



1. INTRODUCTION 

In [8] , the following coupled partial differential equation system was proposed to 
study the invasion by precursor and differentiated cells: 

(1-1) <^ 

[ Vt = Pu{l-f-^), 

where u{x,t) denotes the population densities of the precursor cells. The constant 
d > is the diffusion rate of the cell u which has proliferation rate a > 0, and ki > 
is the carrying capacity of u. The parameter v measures the relative contribution 
that the differentiated cell with population density v(x,t) makes to the carrying 
capacity ki . The cell population density v is limited by its carrying capacity k2 and 
has a max. differentation rate /3 > 0. The model assumes that the differentiated 
cells do not have mobility. 
Letting (see [8]) 

u ^ V ^ ^ j~a 

u — —, V = t = at, X — \ — X 
ki ki y D 

and dropping the hat notation for convenience, system (11.11) is changed into 

Ut = Uxx + u{\ — u — vv) — Au(l — Kv), 



(1.2) 



Vt = Am(1 — Kv), 



where A = ^ and K ^ j^. 

System (jl.ll) or (|1.2p belongs to reaction diffusion systems of degenerate type, 
and such systems have attracted much attention in the research such as epidemics 
and wound healing [QjUniE]. However, system (|1.2|) differs from the above systems 
in the appearance of degenerate reaction terms. In fact, u = coupling with any 
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V — constant consist of a constant solution of (|1.2p . This resembles the combustion 
wave equation considered in [5], however our method in proving the existence of 
the fronts of (|1.2p differs from theirs. 
If the parameters satisfy 

(1.3) 0<v <K, 

then system (|1.2p admits an additional equilibrium: B : (1 ^ "^i ■^) representing 
the state that the spatial domain is successfully invaded. We also separate the 
equilibrium A : (0, 0) from the rest of the line of equilibria, u — 0. The unstable 
equilibrium (0, 0) represents the state before the invasion. 

We are interested in the existence of the wave solutions connecting A with B 
as time and space evolve from —oo to +oo. Setting x + ct,xEM., tE IR+, a 
traveling wave solution to (II. 2p solves 



(1.4) 



— cu^ + u(l — u — vv) — Xu{l — Kv) = 0, 



-cv^ + Xu{l - Kv) = 0, 
with boundary conditions: 

(1.5) (u,i;)(-oo) = (0,0), {u,v){+^) = {l-^,^). 

For the notional convenience we further set 

u 



^ K 



and dropping the bar on u to have 

{u^^ ~ cu^ + u[l - A - (1 - ^)u + {XK - iy)v] = 0, 
-cv^ + X{1 - ^)uil - Kv) = 0, 

and 



(1.7) {u,v){-^)^iO,0), {u,v){+^) = {!,-). 

Numerical investigations '8' strongly suggest that system (|1.8p and (|1.7I) admit 
traveling wave solutions for v = and v = I. When the differentiated cell density 
does not affect the proliferation of the the percursor cells, we have = 0; and when 
the the total cell population contributes to the proliferation carrying capacity, we 
have V — 1. Numerically, however when u = 1, (jl.6p may have non- monotone 
traveling wave solutions and requires a different treatment. Hence in this paper 
we only study the wave solutions for i/ = 0. The system (|1.6p in this case can be 
further reduced to 



(1.8) 



— cu^ + w(l — A — M + XKv) = 0, 
-cv^ + Xu{l - Kv) = 0. 



The computations in [8 shows that the wave may exist for c > 2%/l — A, but 
a rigorous existence proof is stih lacking. We will confirm this observation by a 
mathematical analysis of the model. The system is of cooperative type, we can 
use the monotone iteration scheme developed in 4 for the existence proof. Such 
method reduces the existence of the wave solutions to that of the ordered upper 
and lower solution pairs for ()1.8|) and (|1.7|) . The upper and lower solutions in this 
paper come straightly from two KPP type equations, which are so constructed that 
they have the same decay rate at —oo. Such information is also relevent to the 
monotonicity and uniqueness of the wave solutions. Indeed, since we have a good 
understanding of the decay properties of the solutions at infinities, we then can 
study the properties of the solutions on finite domain, in which the powerful sliding 
domain method (see [IJ) can be used to have the desired results. We remark that 
the methods we used in the proofs of the monotonicity and the uniqueness have 
subtle difference from the ones used in 



2. THE MAIN RESULT. 

In this section we will use monotone iteration method to set up the upper and 
lower solutions for system (II. 8p and (II. 7p . 

Definition 1. A C2(R) x C^{R) function {u{0,viOV, ^ G M is an upper solution 
of dLl]) and dllll) if it satisfies 

— CM^ + u{l — A — M + XKv) < 0, 

-cv^ + Xu{l - Kv) < 
and the boundary conditions 

(2.2) I I (-00) > I ° I ' I I (+^) > 



K 



We can similarly define the lower solution (u, w)(^), ^ G K. by reversing the 
inequalities (|2TT|) and ((Z!2)) . 



The following known result ( iZ] ) is needed in the construction of the upper and 
lower solutions: 

Consider the following form of the KPP equation: 



(2.3) 



u" - CLO' + fiuj)^ 0, 
uj(^oo) — 0, uj{+oo) 



where / G C2([0, 6]) and / > on the open interval (0,6) with /(O) = f{b) = 0, 
/'(O) = a > and f'{b) = -bi < 0. 

Lemma 2. Corresponding to every c > 2\/E, system has a unique (up to a 

translation of the origin) monotonically increasing traveling wave solution w(^) for 
^ G M. The traveling wave solution uj has the following asymptotic behaviors: 
For the wave solution with non-critical speed c > 2\/d,, we have 



(2.4) uj{^) = Qu^e" ^"^ + o(e^^^T^«) as ^ 



(2.5) uj{^)^b~b^e ^ « + o(e 2 as ^ ^ +00, 

where and are positive constants. 

For the wave with critical speed c — 1\f^, we have 

(2.6) = 4Ce^« + o(Ce^«) asi^ -00, 

(2.7) uj[i) = b - b,e^^-'^^^^ + o(e(^-^'^)«) as ^ ^ +00, 
where the constant dc is negative, be is positive. 

We next consider the following version of the KPP system: 



(2.8) 



uj" - cw' + (1 - A)w(l - w) = 0, 
w(— 00) = 0, aj(+oo) = 1. 



According to Lemma [21 for every c > 2v]— A system (|2.8p has a unique (up to 
a translation of the origin) monotone solution u{£,), ^ € M. Now fix this u{£,) and 
consider the equation 

(2.9) - cv^ + Au(l ~Kv)^0 



For each fixed c > 2Vl — A and the corresponding u{£,), (|2.9p has a solution 

(2.10) ii(e) = ^(l-e-^^-^--"(^)''^). 
We next compare £t(^) with v{£,) for ^ g K. 

Lemma 3. There exists a (^i > such that if 

(2.11) 0<A< 



2 + K{l + V2)' 
we have 

Proof. According to Lemma [21 the wave solution {t(^) to (|2.8p has the following 
asymptotic behaviors: 
For c > 2yi~A, 



(2.12) ?i(^) = a^^e ^ « + o(e s «) as ^ ^ -00, 



(2.13) w(0 = & - fo^e 5 « + o(e ^ «) as ^ ^ +00, 



and Uuj, bui are positive constants. 
For c = 2v'l — A, we have 



(2.14) u{0 = dc^e"^^ + o{£,e^^^) as ^ -oo, 

(2.15) u{0 = b - 6ee(i-^)^« + o(e(i-^)^«) as ^ ^ +cx). 

where the constant dc is negative, be is positive. 

We now study the asymptotics of the function v{^). Formulas (|2.12l) and (I2.14p 
imply that 

/ u{s)ds as^— >— oo. 

We can then expand 

(2.16) e-^^-^-~ = 1--K u{s)ds + o(( u{s)dsy). 

^ J — oo J — oo 



A further expanding of (|2.16p for ^ — > — oo and for c > 2^/1 — A 
(2.17) 

" (c- Vc2-4(1-A))c 



and for c = 2Vl — A, 

(2.18) 1 - e--^^'-^ ^^^^"^ = ^^^(4^ - T^)^^' + o(Ce^«). 
As for ^ > sufficiently large we have 

uls)ds 

(2.19) lim / = lim u{£,) = 1, 

5— S- + 00 ^ {— >+oo 

therefore 



(2.20) viO ^ K{1~ e--^^) +oie~-^^) as^^+oo. 

We next show 



(2.21) _X^^c_-^lTWER_ 



or equivalently 

2c(l - A) 



:+ Vc2+4(1-A) 



> XK. 



Setting q{c) = .^''^ , then it is easy to see that q(c) increases as c 

does. Hence 

2(1 -A) 



5(c) > 5(2%/r^) = 



We therefore require < A < 2+k(i+\/2) ^^'^^ (I2.2ip . 

We now shift u(^). Since (|2.8p is shifting invariant, + (),£, ^ R is also a 
soKition for any C e M. It then foUows from (|2.12p for c > \/l — A, 

c-y/c^ -4(1-A) ^ c- v^c^ -4(1-A) ^ c-\/c^-4(l-A) ^ 

(2.22) it(e + C) = fl^e^^T ^e^^T « + o(e^^T «) as C ^ -oo; 



and for c = Vl — A, 

(2.23) u{£, + C) = 4(^ + Oe^^'^e'^^i + o{£,e^^^) as ^ ^ -oo. 



c—-,/c2— 4(1 — A) 

If we choose C > sufficiently large, the positiveness of — - — ^ C, and 



y/1 — AC implies that if c = 2\/l — A, 



2vT^' 1-A' 
and if c > 2^/1 - A, 



c-v'c^-4(i-A) ^ 2Aa^ 



(c- v/c2-4(l-A))c 

It then follows from ((2?2T1) . (|Z;22)) and ([S^ that there exists a iV > sufficiently 
large, 

^u{£_ + C) > v{£,) for £_ e (-oo, -iV] U [N, +oo), 

and for ^ e [— A^, A^], since m and v are both monotonically increasing on M we can 
further shift + C) to the left at most 27V units to have ^u(£, + C) > ^(0, C ^ K. 
Hence there exists a finite Ci ^ such that the conclusion of the Lemma holds. □ 

Now we write u{0 = + Co), ^ G M and let v{^) be defined in ([2J0|) . We 
remark here that the computation of v(^) still uses u(^). 

Lemma 4. Assume t/ie conditions in Lemma\^then (t2,'u)(^), f S M defines an 
upper solution for \1.8]) and |i. 7] ). 

Proof. We can easily verify that (m, u)(^) satisfies the boundary conditions 
For the u component we have 

u" - cu' + u(l - A - M + \Kv) 

= u[l-X-u + \Kv-{l~\){l~u)] 

= -XKu{^u-v) <0 

The last inequality follows from the previous Lemma. 
As for the v component, for each u, we have 



-cv^ + Am(1 - Kv) = 0. 



□ 



We next set up the lower solution for (jl.8|) and (|1.7p . 

For a fixed I > we consider another version of the KPP system: 



(2.24) 



w" - cw' + (1 - X)w{l - ^i-r^w) = 0, 



w(-oo) = 0, u;(+oo) = j^j^i < 1- 



Then for any c > 2Vl — A, (|2.24p has correspondingly a unique wave solution u(^), 
We define 

The next Lemma gives the relation between and C G H^- 
Lemma 5. There exists a Ci > such that 

(2.25) ^H^-Ci)<v{0 CeK. 

Proof. The proof is similar to that of Lemma [3l Noting as ^ — > +oo, u{^) — 
"^T^^ITI ^ Hence we do not need condition p. lip here. □ 

We denote u{£,) = u(C - Ci), C e K. Then 
Lemma 6. Such defined (w, w)(0; C G consists of a lower solution for il.8\) and 

Proof. One the boundary we have 

{u,v){^^) = (0,0), iu,v){+^) = (y^^' J^) < (1' J^)- 
and for the u component, 

u" - cm' + m(1 - A - u + XKv) 
= u" - cu' + (1 - A)u(l - ij^w) - (1 - A)m(1 - ^t^m) 
+u(l - A - M + Aifw) 

= m[1-A-u + AXi;-(1-A) + (1-A + 0u] >0 

due to the last Lemma. 

Noting that v solves the equation 

-cv^ + Ait(l - Kv) = 0, 

it satsfies the inequality trivially. □ 

Lemma 7. The upper and lower solutions are ordered 

(2.26) {u,v)iO>iu,Em, eeK. 



Proof. For each fixed c > 2\/l — A , if u(0 soives tlie system (I2.24[) tlien tlie function 
HO = solves ^jSI- Hence it follows that u{C) > HO hence u{C) > u{0 

for aU e e R. 

By the definition of v{^) and v{Oi '^^ have 



> 

Hence the conclusion of the Lemma holds. 



□ 



Theorem 8. Let the parameters satisfy i2.11\) . then for each c > 2\/l — X, system 
lil.8\) and |_?. 7| ) has a unique (up to a translation) strictly monotonically increasing 
traveling wave solution, while for < c < 2\/l — X, there is no monotonic traveling 
wave. The Traveling wave solution has the following asymptotic behaviors: 
For c = 2Vl - A 



(2.27) 



(0 = 




+ o(e^i-^«) 




j_ 

K 



C2ie 



C22e 




/ o(e 



and for c > 2\/l — A 
(2.29) 



where 011,012,021,022,(^21,^22 > 0, dii,di2 < 0. 



as e — !> —00, 





as e 



-00; 



as e —00, 



as e +c 



Proof. Noting that between the upper and lower solutions, there is no equilibrium 
other than (0, 0) and (1, -^) of system p.Sp and (|1.7p . Hence the monotone iteration 
scheme developed in is still applicable. Such monotone iteration scheme reduces 
the existence of the traveling wave solutions to that of the ordered upper and 
lower solution pairs, the existence of the traveling waves then follows by Lemma 
[6l Lemma |4] and Lemma [3 and by [4] , such obtained traveling wave solutions 
are nondecreasing. While for c < 2-\/l — A it is easy to verify, by analyzing the 
equilibrium (0,0) that the nontrivial bounded solutions of (jl.Sp are oscillatory. 



We next show that the wave solutions are strictly monotonically increasing on 



For any fixed c > Vl — A, let {uc,Vc){£,) be the corresponding traveling wave 
solution, and (u'i(^), u)2(C)) be its derivative. Then (u'i(^), w;2(C)) > for ^ € R, 
and (wi(^), ?«2(C)) satisfies the following systems 



(2.31) 



wi^^^ — cwi^^ + (1 — A — 2uc + XKvc)wi + \KucW2 ~ 0, 

— CW2,5 + A(l — KVc)wi — \KUcW2 — 0, 



(wi, W2)(±oo) = 0. 
It then follows that 



(2.32) 



wi^jj — cwi_5 + (1 — A — 2uc + XKvc)wi < 0, 

— CW2.^ + A(l — KVc)wi — \KUcW2 = 0, 



(wi,W2)(±00) = 0. 

Applying the Maximum Principle to the first inequality of ()2.32p . we immediately 
conclude that 'U'i(^) > for ^ G K. Thus Uc(0 is strictly monotonically increasing. 

The strict monotonicity of fc(C) comes from (|1.8p . Since Uc{C) > for all ^ € M, 
and for such Uc(0 have 



K 



(1 



then it follows that W2{C) = "v'dO > 0, ^ G M. This shows that the wave solution 
{uctVc) is strictly monotonically increasing. 

We then derive the asymptotics of the wave solutions at ±oo. Noting that the 
upper and lower solutions have the same exponential decay rate at — oo, (|2.27p and 
(|2.29l) come directly from comparison. 

We next study the asymptotics of the function {wi,W2)[C) at +oo, recalling 
that {wi,W2)[C) — {uc,Vc)'{C) and satisfies the system ()2.3H) . Since this system 
is hyperbolic at +oo, {wi,W2) approaches (0,0) exponentially. We will derive the 
exact exponential rate. 

The limit equation at +oo of system (|2.3ip is 



cw'^ t — Wi + \Kw2 = 0, 



— J — \KucW2 = 
Since the second equation is decoupled from the system, we immediately have 



w. 



Plugging the above into the first equation yields a bounded solution (up to the first 
order approximation) of the form 



v+{0^A,e-^i^A2 



(t) = Aie <= ^ + A2e 2 
By roughness of exponential dichotomy [?], we have 



MO ) V ^e^« / V 
where n is either — ^ or ^^^^^S, 

Integrating the above from to +00, and comparing the decay rates of (uc, Vc)i£,) 
with that of the upper solution {u, we have (|2.28l) and (|2.30p . 

On the uniqueness of the travehng wave solution for every c > 2\/l — A, we 
only prove the conclusion for traveling wave solutions with asymptotic rates given 
in (|2.29p and (|2.30l) since the other case can be proved similarly. Let Ui{^) = 
(wi,wi)(^) and t/2(0 = ('^2,W2)(^) be two traveling wave solutions of system (ITT 



and (|1.7p with the same speed c > 2y/l — A. There exist positive constants Aij, 
Bij, i,j — 1,2 and a large number > such that for ^ < — A^, 



c + Ve^-4(l->) ^ c + v'c^-4(l-A) , 



(2.33) t/i(0=| |e ^ ^ + o(e ^ ^) 

^12 

A21 

A22 



(2.34) C/2(C)= e ^ ^ ^ + o{e ^ ' 



and for f > TV, 

(2.35) f/i(0 = 

(2.36) ;72(0 = 



The traveling wave solutions of system (|1.8l) - (|1.7p are translation invariant, thus 
for any 6* > 0, f/f (^) := C/i(f + 9) is also a traveling wave solution of pTS jl -pTT l) . 
By (|2.33p and (|2.35p . the solution Ui{^ + 6) has the asymptotics 




(2.37) C/f (0 = 

for ^ < -iV; 



/ o+\/c2-4(l-A) , 
. c+Vc^-4(l-A) , 

V Aiae 2 ^ 



c + Vc2-4(l-A) ^ o+\/c2-4(l-A) , 

5 « + o(e 5 « 



(2.38) C/f(e) i ^ , . , 

A / ^ ' 

for C > A^- 

Choosing 6* > large enough such that 



c + v'c2-4(l-A) a 

(2.39) Aiie ^ " > ^21, 



(2.40) 



(2.41) 



(2.42) Si2e--^^ < B22. 
then one has for ^ e (—00, —N\j[N, +00), 

(2.43) UliO > C/2(0- 

We now consider system (|1.8p on [— iV, +A^]. There are two possibihties: 

Case 1. Suppose we already ahve C/f (C) > f^2(0 O'^ [~-^: then the function 

M/(^) = (u;i(^),u;2(0)'^ := C^f(0 " ^^2(0 > and satisfies for some G (0,1), 

* = 1,2, 



(2.44) 








H/(-iV)>0, M/(+7V)>0. ^e{~N,N), 
where the matrix M is given by 



(2.45) 

M{wi,W2) = 



1 - A - 2(m2 + ClWl) + \K{V2 + C2W2), \K{U2 + Cl""^!) 



A(l-i^(«2+C2W2)), 



'\K{U2 + C2W1 



Since u;i(^) > 0, C G [-N,N] and \K{u2 + Ci^i) > 0, then we have on ^ e 



(2.46) 



cw[ + [1 - A - 2(m2 + ClW'l) + \K{V2 + (2^2)]^! + \K{U2 + C2Wi)w2 = 0, 



u;i(-iV)>0, wi(iV)>0. 

The Maximum Principle then implies that wi(^) > on [—N, N]. We then move 
to the second equation of (|2.44p . We have 

(2.47) 

-cw'^- XK{u2 + C2Wi)w2 = -Xil-K{v2 + C2W2))wi<0, ^e[~N,N], 

W2{-N) > 0, W2{N) > 0. 

The strict inequality comes from the fact that 1^2 (C) ^ i'2(C) + C2W'2(0 1^ < ^ 
for ^ e [-N,N]. It then follows that W2{$.) > for ^ G [-N,N]. For if there is a 
^ e {—N,N) such that W2(C) = O.Then u'2 takes local minimum at then the left 
hand side of the first inequality of p.47p is zero at We then have a contradiction. 

Case two. We may suppose that there is some point in {—N, N) such that one 
of the components, say the j-th component, satisfies (J7i(C))j < {U2{£,))j at that 



point, j ~ 1 or 2. We then increase 9, that shifts C/i(C) further left, so that 
Uf{-N) > U2{-N), Uf{N) > U2{N). By the monotonicity of C/f and C/2, we can 
find a 6* e (0,27V) such that in the interval (-iV,iV)^we have U%£, + 9) > U2{£.)- 
Shifting ?7f (^+^) back until one component of C/f first touches its counterpart 

of U2{£,) at some point ^ e [—N, N]. We then return back to case 1 again, where it 
has been shown that this is impossible. Therefore, we must have 

uliO > U2{0 

for all ^ e M, where 9 is the one chosen by means of (|2.39l) - (|2.42p as described 
above. 

Now, decrease 9 until one of the following situations happens. 

1. There exists a ^ > 0, such that Uf{£,) = U2{(.). In this case we have finished 
the proof. 

2. There exists a, 9 > and ^1 e R, such that one of the components of and 
U2 are equal there; and for all ^ e M, we have C/f (^) > U2{^). On applying the 
Maximum Principle on R and use the same argument as we did for case 1. We see 
this is impossible. 

Consequently, in either situation, there exists a 9 > 0, such that 

uhO ^ U2{0- 

for aU C e R. □ 
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